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! ABSTRACT 

I 

Axially-dependent  Vlasov-Poisson  equilibria  are  studied 
for  monoenergetic  unneutralized  electron  beams  injected 


along  magnetic  field  lines  into  a vacuum  half-space  and 

j 

partially  or  totally  reflected  by  the  beam  space  charge. 


For  simplicity,  the  entering  beam  electrons  are  assumed 
to  be  nonrelativistic  and  to  have  a rectangular  distribution 
in  momentum  pitch  angles  with  respect  to  the  magnetic 
field.  Electron  density  at  any  point  along  a field  line 
is  first  expressed  in  terms  of  a potential  for  parallel 
motion  (assuming  conservation  of  magnetic  moment) , and  this  is 
then  substituted  into  Poisson's  equation  to  obtain  a nonlinear 
partial  differential  equation  for  the  space  charge  potential. 

The  one-dimensional  problem,  in  which  only  axial  electric 
fields  are  considered,  is  reduced  to  quadrature  and  an  exact 
solution  is  obtained  in  the  singular  case  where  the  beam 
has  no  spread  in  pitch  angles.  The  two-dimensional  problem, 
in  which  radial  electric  fields  transfer  energy  from  the 
axial  streaming  motion  to  azimuthal  E*B  rotation  to  assist 
the  axial  reflection,  has  been  reduced  to  a nonlinear 
partial  differential  suitable  for  computational  solution. 
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I . INTRODUCTION 

Virtual  cathode  formation  in  vacuum-propagating  beams'^"  may  be 

more  widespread  than  previously  recognized.  The  "spin  death"  of 

parallel  beam  energy  by  Erxl*z  rotation  in  front  of  a density  pileup, 

may  assist  or  steepen  the  axial  density  gradient,  since  the  rotation 

speed  tends  toward  Brillouin  flow,  a)=u)  / 2.  It  may  be  that  the 

c 

2 2 

propagation  limit  (on  the  beam  density)  is  not  w 'v-w  in  the  applied 

pc 

2 

or  total  magnetic  field,  but  rather  to  (in  the  presence  of  self- 

p max 

field  pileup)  small  enough  that  cE^/B^cVq  (Vg=beam  electron  speed) , 

with  given  self-consistently  from  the  charge  distribution. 

In  what  follows,  self-consistent  BGK-like  axially  dependent 
2 

equilibria  will  be  described,  in  which  the  space  charge  effects 

determine  the  axial  extent  and  profile  of  the  beam.  Generalization 

to  axially  nonuniform  magnetic  fields  is  relatively  straightforward. 

Sections  III  and  IV  treat  the  problem  in  a purely  one-dimensional 

way,  keeping  the  electrostatic  effects  of  Ez  electric  fields  but  not 

the  effects  of  beam  rotation  caused  by  Er.  Section  V generalizes 

these  results  to  include  E . Section  VI  summarizes  the  conclusions 

r 

of  this  introductory  work. 


II.  APPROXIMATIONS  AND  ASSUMPTIONS 


To  calculate  the  electron  density  n(s)  at  a point  s along 
magnetic  field  lines  in  terms  of  an  effective  potential  U(s)  for  the 
parallel  motion,  it  is  assumed  that  all  changes  of  U(s)  and  n(s)  with 
time  occur  much  slower  than  the  time  required  for  a typical  particle 
transit  from  s=0  to  the  end  of  the  system  and  back,  and  that  the  beam 
is  on,  supplying  electrons  at  a constant  rate,  for  a time  much  longer  than 
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the  transit  time.  It  is  also  assumed  that  the  beam  source  absorbs  all 
reflected  electrons,  that  is,  reflexing  is  neglected.  This  "fast  bounce 
approximation"  then  allows  use  of  the  steady  state  Vlasov  equation,  in 
its  integrated  form,  to  calculate  n(s)  in  terms  of  the  momentum 
distribution  f (p,s=0)  at  s=0  and  the  generalized  potential  U(s).^ 

The  potential  for  electron  parallel  motion  is  taken  to  have  the 
form  U(s)=yB(s)-e4»(s) , where  y is  the  magnetic  moment  invariant, 

B(s)  is  the  magnetic  field  strength,  and  <t>( s)  includes  all  effects  due 
to  the  electric  field.  $(s)  may  or  may  not  be  approximated  as  a 
purely  electrostatic  potential.  The  electrostatic  potential  <f>(s) 
defined  by  E)(  =-9<J>/3s  enters  into  $(s),  but  e$(s)  in  general  also 
contains  the  energy  in  ExB  drift  motion  (taken  out  of  parallel  motion.) 

(In  general  the  variables  U,  B,  n,  <J>,  etc.,  may  also  have  variation 
with  coordinates  normal  to  field  lines,  provided  the  scale  lengths 
for  such  dependence  are  larger  than  a gyroradius.) 

It  will  be  assumed  that  there  is  some  point,  s=M,  such  that  for  s>M, 
no  particles  are  reflected  at  s.  This  point  may  be  at  infinity,  but 
it  is  convenient  to  keep  it  in  the  analysis. 

For  simplicity  in  evaluating  the  integrals  involving  f(p,s*0), 
the  electron  beam  distribution  in  momentum  will  be  assumed  monoenergetic 
with  a rectangular  distribution  in  pitch  angles.  (More  general  angle 
distributions  are  also  tractable.)  Also  for  simplicity,  we  take  the 
case  where  9B/9s=0,  i.e.,  where  the  reflection  of  beam  electrons  is 
entirely  due  to  electric  field  effects  and  not  to  magnetic  mirroring. 

(This  is  easily  generalized  for  nonrelativist ic  motion,  and  changes 
only  the  boundary  of  integration  of  the  integral  expression  for  n(s).) 

Finally,  since  we  want  only  to  show  the  basic  physics  of  beam 
self-reflection,  we  consider  nonrelativistic  beam  momenta.  Generalization 
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to  relativistic  motion  is  straightforward  but  mixes  p)(  and  pA  through 
the  factor 

J,  . . 2 . 2,  . 2 2 

Y = T1  + (p„  + Pj.)/m0c  , 

which  occurs  in  the  n(s)  integrals. 

The  collisionless,  magnetic-moment-preserving  behavior  assumed  here 

implies  that  instabilities  do  not  disrupt  the  electron  motion.  This 

requirement  may  be  satisfied  since  the  virtual  cathode  formation  will 

cause  inhomogeneity  of  the  density  on  the  scale  of  one  beam  Debye  length. 

Thus,  even  though  reflection  of  a beam  with  small  parallel  momentum 

spread  generally  allows  counterstreaming  instability  in  homogeneous 

plasmas,  this  instability  may  not  be  able  to  grow  to  significant 

amplitude  in  the  space  provided,  all  the  more  because  density 

gradients  absorb  the  unstable  waves.  Radial  motion  due  to  E xB 

z u 

is  neglected.  This  additional  self-pinching  due  to  E must  be  offset 

z 

by  centrifugal  effects  arising  from  the  curvature  of  the  helical  magnetic 
field  lines. 

III.  DENSITY  IN  TERMS  OF  POTENTIAL 

The  density  at  s is  comprised  of  forward-going  and  reflected 
electrons: 

n(s)  = n+(s)+n_(s)  . (1) 

I 

The  velocity  distribution  of  the  forward-going  portion  is  specified 
at  s»0: 

f+(v„  .vJL,s-0)  = N<5(v-  vq)H(v)H  (y0-v±)  (2) 

where  H(x)  is  the  Heaviside  function:  H(x)*l  for  x>0,  -0  for  x<0,  and 


Yq  is  a maximum  value  of  vf.  This  distribution  is  nonzero  on  an  arc 
of  a circle  in  v-space  at  s=0.  See  Fig.  1. 
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The  velocity  distribution  of  the  reflected  electrons,  f (v|(  ,vA,0) 

follows  from  reflection  conditions  and  so  is  not  specified  a priori. 

4 

One  can  now  use  the  integrated  moment  form  of  the  one-dimensional 
steady  state  Vlasov  equation  to  express  n(s)  in  terms  of  f+(v,0)  and 


4>(s)  : 


n(s)  = 2ir[B(s)/B(0)  ] 


>oo 

ydy  xdx(x2-x2)  ^ 2f(x,y,s=0) 

Jo  Jj  s 


* 5 v(|  , y = vx  , 


where 


is  the  square  of  the  'escape'  velocity  required  to  reach  point  s if  x is 
monotonically  increasing  in  s,and  the  integration  region  £ in  x will  be 
discussed  presently.  Let 

& =Sup  x , (5) 

S s'<s  S 

be  the  maximum  value  taken  on  by  x , at  points  between  s'=0  and  s'=s,  i.e., 

s 

the  "escape"  velocity  required  to  reach  point  s even  if  xg  is  not  monotonic  in  s. 

r\j 

Particles  with  parallel  velocity  x>xg  at  s=0  contribute 

% 

to  the  density  at  s;  those  with  0<x<xg  have  been  reflected 

at  some  s'<s  and  do  not  appear  at  s.  Because  of  the  beamlike  nature 

of  f+(v,0),  there  are  no  electrons  with  positive  x less  than  some  xQ  (i.e., 

with  y greater  than  some  y^) ; thus  for  s less  than  some  Sq  there  are 

no  turning  points.  For  s>Sq  there  are  turning  points  until  s>M.  The 

point  Sq  is  specified  by  its  potential,  i.e., 

2 2 

x = x 

s0  0 


with  x given  by  Eq.  (4)  with  s replaced  by  s_.  The  minimum  parallel 
s0  u 

velocity,  Xq,  at  the  s=0  reference  plane,  is  a parameter  of  the  distribution 


function.  There  are  no  turning  points  at  points  s<Sq. 

For  points  s>Sq,  the  lower  boundary  of  the  x- integration  in 

% 

Eq.  (3)  is  xg  as  given  in  Eq.  (5).  For  s<Sq,  the  minimum  is  simply 
Xq.  In  both  cases  the  maximum  can  be  taken  as  Vq,  as  is  clear  from 
Fig.  1.  Thus,  when  B(s)=B(0),  one  has  for  the  forward-going  electrons 


wj  t 2 2.  2 

Minys’y0  rv0 


n+(s)  = j N j d (y2)  j 


' ° d(x2)  /T~l 

— <$(/x  +y  -v_) 

^2  2 fl  2 u 

_x  or  x„  vx  -x 
s 0 s 


fl  2 

= 27rv„N/v  -x  for  s>s„  , 
0 0s  0 


= 2Trv0N[4^-7-  X2-y2-x2]  for  s<sQ 


2_  2 2 


where  ys=vQ~xs*  In  2,  the  forward-going  density  n+(s)  is  shown 

2 

vs.  xg,  i.e.,  vs.  the  normalized  potential  difference 


ij>=(2e/mv0)  [4>(0)-$(s)  ] 


The  maximum  value  of  n+  in  Fig.  2 occurs  at  ^o=x(r  v0‘ 

The  density  of  reflected  particles  is  calculated  similarly. 

One  assumes  there  is  an  s=M  such  that  no  reflections  occur  beyond  M. 
Then  n_(s)  is  comprised  of  particles  that  go  past  s but  not  past  M: 

V>0  2 

- f » d(y2)  J_  a(^V-Vn)  o 

2 j 7 n 7 0 


J2  >2/22 
y_  x /x  -x 

m s s 

= 2irv  n/x2-x2  for  s„<s<M 
0ms  0 


Note  that 


2itv  N [ ^-xf-/x2-x2  ] for  s<sn 
U m s 0 s 0 


x2-x2  = (2e/m)  [ 4> ( M)  — 4> ( s ) ] . 
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2 2 2 2 2 2 

For  convenience,  let  ili=x  /v.,  il»  =x  /v-.  The  physical 

s 0 0 0 0 m m 0 

range  of  the  "independent  variable"  ip  is  then  from  zero  ( if; (0)  =0)  to 

ik.  (M)  =d»  ).  The  density  at  s is  then 
M m 

n(s)=n,  (s)+n  (s)=2irNv^ [/l-ip-S\pn-\p  + -\p-i/\p  — ] (14) 

+ - U U m u 

with  any  radical  dropped  once  it  becomes  imaginary.  The  density  n(s) 
is  shown  vs.  ip  in  Fig.  3. 

IV.  SELF  CONSISTENT  SOLUTION  FOR  AB=0 


The  equation  for  n(s)  in  terms  of  \p(s ) is  now  coupled  to  a 
Poisson-like  equation  for  ip(s)  whose  right-hand  side  is  n(s) . 
The  ^suiting  equation  has  the  form 


/ (4>)  = dip/ dr) 


(15) 


where  u>  is  a functional  of  (and  in  general,  its  derivatives) 

P 

and  f is  the  Laplacian  operator.  The  simplest  model  is  a purely 

one-dimensional  one  with  4>(s)=4> (s)  (i.e.,  purely  electrostatic  potential) 

2 2 1/2 

and  /=d  /ds  . In  terms  of  the  scaled  distance  z=s*(2  v/v  the 

equation  is 


[l-2/ip7+Jip~  ] 
0 m 


0<tp<\p 


(16) 


2 2 

where  ^qHxq^vq  an<*  where  the  normalization  N has  been  evaluated  in  terms  of 

- ...  2 2 
the  density  n^  at  s=0  and  subsumed  into  the  definition  of  z via  cd^Q=4irnQe  /m. 

It  is  not  known  whether  this  nonlinear  equation  always  has  unique 

solutions.  When  the  spread  in  parallel  velocities  of  the  beam  is  small 

we  can  cast  this  equation  in  a soluble  form  if  all  electrons  are  reflected 

1 2 

somewhere:  then  \p  =1,  i.e.,  e [ <J>  (0)  — (M)  ] = -r-  mv. , and  the  right-hand  side  of 
m L u 


(16)  is 


9 


■/i— \|i — iApQ— xp 

i-^ 

For  small  angle  spread  we  may  write  i|/g=l-e  (see  Fig.  1)  where 

2 

in  terms  of  the  pitch  angle  spread  eQ,  e=eQ/2.  Expanding  in  e, 
the  right-hand  side  becomes  approximately 


/l-i p 


for  1-\|j>>£,  where  a-»-l  as  e-K>. 
The  resulting  equation, 


,2. 

d 1 1)  _ _ q 

O 

dz  /l-iji 


(for  l-i p»e) 


(17) 


is  a turning  point  problem  and  can  be  solved  by  quadrature: 


z = [d-v^)(i+2v/iq;)1/2] 


(18) 


giving  ( z)  implicitly  as  z ( \p)  . In  particular,  the  condition  that 
no  particles  are  reflected  beyond  some  point  s=M  (and  that  all  are 
reflected  somewhere)  gives  i|)  or  M,  subject  to  'i'm$l: 

M = [(l-^T)(l+2/iT)1/2l  • (19) 

3a  m m 

See  Fig.  4,  for  which  the  limiting  value  a=l  has  been  used.  In  this 
case,  the  infinite  charge  density  at  s=M  leads  to  the  discontinuity 
in  the  slope  of  f at  M in  Fig.  4.  When  e^O  this  discontinuity  disappears 
and  i p has  a simple  maximum  in  s. 

In  the  limit  of  no  angle  spread,  all  the  electrons  are  reflected 
at  this  same  point.  One  can  then  use  Eqs.  (14)  and  (18)  to  obtain  n(s) : 


n(s)  = l/u(s)  u(s)  = Jl-ip 

u3  - | u2  (l-9z2)  = 0 


(20) 


(21) 


V.  EFFECT  OF  RADIAL  ELECTROSTATIC  FIELD 


Although  the  treatment  in  Eqs.  (16)-(25)  is  strictly  one-dimensional, 
i.e.,  involves  only  electric  fields  along  s,  a formal  generalization 
is  also  possible  when  electron  parallel  streaming  energy  is  converted 
to  azimuthal  EQxB  rotation  due  to  the  radial  electrostatic  self-field, 
EgCr.s).  Since  any  increase  in  the  energy  of  this  azimuthal  motion, 

j m(cE0/B)2  , 

m.’st  come  from  the  initial  streaming  energy,  it  can  be  treated  as 
part  of  the  effective  potential  for  parallel  motion,  U(s) , which  now 
becomes 


U(r,s)  =|  y2B01B(r,s)-e<j)(r,s)  +f  c2  [EQ(r , s) /B(r , s)  ]2 


(26) 


As  mentioned  in  Sec.  I,  this  can  be  written 


U(r,s)  y2BQ1B(r,s)-e<|>(r,s) 


(27) 


where  $ now  includes  all  effects  due  to  the  electric  field,  including 
the  radial  field,  and  where  the  magnetic  moment, 


is  assumed  conserved. 

For  arbitrarily  large  beam  density,  azimuthal  current  now  makes 


B(r,s)  nonuniform: 


Vx^(r.s)  n(r,s) 

c 


cE(r,s) 
B(r , s) 


(28) 


and  the  natural  coordinates  ( n i S)  become  noncylindrical,  complicating 
2 

the  V operator.  However,  to  demonstrate  the  physics  of  the  electrostatic 

effects  alone,  we  may  now  restrict  the  remaining  discussion  in  this 

section  to  the  case  of  small  beam  radius,  low  density,  or  strong  I^q. 

Writing  Eq.  (28)  in  dimensional  form  with  S(r,s)=§ n+§  (r,s),  one  can 

u s 


show  that 


2 


-1 
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with  r,  the  beam  radius  and  w n the  plasma  frequency  based  on  the  beam 
b pu 

density  at  s=0.  Thus  in  Eqs.  (30)-(35),  we  shall  require 

2 

r,  a) 

_E 

c w 


« 1 


(30) 


so  that  the  magnetic  self-fields  are  weak,  while  the  electrostatic 
self-fields  are  still  strong  enough  to  be  important.  In  this  case, 
because  of  the  electrostatic  nature  of  the  problem,  one  has 


EQ(r,s)  = - — <j>(r,s) 


and  straight  magnetic  field  lines.  Thus  we  can  write 

2 

(?(r,s)=<j>(r,s)-(mc2/2eB2)(||:) 


(31) 


(32) 


for  the  case  B(r ,s)=const.=BQ. 

Because  reflection  of  electrons  depends  only  on  the  potential  for 
parallel  motion,  an  azimuthally  symmetric  nonneutral  distribution  of 
electrons  will  have  its  density,  n(r,s)  given  by  Eq.  (3),  as  before, 
except  for  the  modification  in  the  definition  of  <J>.  For  the 
distribution  function  of  Eq.  (2),  the  results  of  Eqs.  (6)-(14) 
remain  valid.  A self-consistent  equilibrium  for  straight  magnetic  field 
lines  can  now  be  obtained  in  the  form  of  Eq.  (15);  letting 


2e  , . 

\p  = — 2 <Hr*£0 
mv0 

be  the  scaled  electrostatic  potential,  and  defining 


(33) 


[Vq/ (2eBg/mc) ] =b  , (34) 

2 

we  may  simply  replace  ^ by  Y+bOty/ar)  in  the  right-hand  side  of  Eq. 

2 

(16).  Note  that  b is  just  1/4  of  the  (gyroradius)  , so  that  the 
natural  radius  scaling  is  of  the  order  of  the  gyroradius.  Eq.  (16) 
then  becomes,  for  the  2-dimensional  electrostatic  problem, 

(rr  )r  h (r  If ) + M = -|^-*-b(3*/3r)2 

\2o)  / 32 

N p0/  / s 7 

-2/cos  eQ— ip— b (3ip /3r) 


axz  [vp-t-b  (dip/d r)  2 ] - ty+b (a^ /a r ) 2 ] 


-1 


cos^Gg-g  +V>laxztiJ)+b(a4//ar) 


(35) 


WM 


where  0q  is  the  initial  spread  in  pitch  angles  of  beam  electrons,  gHb(3t|//3r) 
and  where  the  operation  Max^  is  taken  at  fixed  r.  The  nonlinear  partial 
differential  equation  (35)  appears  suitable  for  moderately  simple 
computer  solution  by  overrelaxation  methods. 

The  condition  (30),  together  with  the  requirement  that  the  electron 
gyroradius  be  small  compared  with  the  beam  size  r^,  give  the  requirement 

v0/c«(mc/u)p)2  » -V 

where  oo  EeB^/mc. 
c U 


VI . CONCLUSIONS 


In  the  one-dimensional  case  of  purely  electrostatic  repulsion 

(Ez>0,  Er=0) , self  reflection  of  the  entire  electron  population  occurs 

in  a finite  distance  even  if  there  is  no  pitch  angle  spread.  For 

small  or  zero  spread,  this  virtual  cathode  occurs  at  a distance  L'vv^/bu) 

'v  0 

i.e.,  at  less  than  one  beam  Debye  length  based  on  the  initial  kinetic 
energy  per  particle.  For  a cold  beam,  all  the  beam  electrons  are 
reflected  at  the  same  point,  namely  jit  this  distance  L.  The  density 
profile  is  then  given  by  n(s)=l/u  with  u the  solution  of  the  cubic 
equation  (21).  The  purely  electrostatic  one-dimensional  problem  with 
angle  spread  has  been  reduced  to  quadrature,  as  given  by  Eqs.  (24)  and 
(25). 

For  the  case  with  E^O,  the  scaled  electrostatic  potential 
2 

i/)=2e(J)/mvQ  is  given  by  the  parabolic  equation  (35).  This  equation  is 
suitable  for  moderately  simple  computer  solution  by  the  method  of 


overrelaxation. 
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FIGURE  CAPTIONS 

Fig.  1 Reference-plane  velocity  space,  x«v(|  > y*v±,  at  spatial  point  s-0. 
Beam  input  velocity  distribution  Is  nonzero  on  the  spherical 
cap  with  x>Xq.  Beam  electrons  with  Xq<x<xs  are  reflected 
somewhere  in  the  interval  (0,s)  and  do  not  contribute  to  n(s). 
Reflected  particles  (x<0)  also  contribute  to  n(s).  There  may 

<x<Vq  are  never  reflected. 
In  this  case,  -Vq<x<-x^  does  not  contribute  to  n(s) . 

rig.  2 Forward-going  electron  density  n+(s)  at  s [Eqs.  (7)  and  (8)), 
vs.  the  potential  difference  [Eq.  (9)]  between  point  s and 
s-0;  shown  for  beam  pitch  angle  spread  0^  at  s=0  given  by 
sln2e0=y2/v2-0,  0.1,  0.3,  0.5. 

Fig.  3 Total  electron  density  n(s)  at  s,  vs.  the  normalized  potential 
difference  ^ between  point  s and  s-0  [Eq.  (14)1.  Shown  for 
ij/Q-0. 7,  i^-0. 95.  The  lower,  dashed  curve  describes  the 
never-reflected  electrons  at  s>M,  i.e.,  beyond  the  potential 
maximum. 

Fig.  4 iHz)  for  z<M,  from  Eq.  (18). 


be  an  x^<Vq  such  that  electrons  with  x^ 
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